Interplay between pairing and exchange in small metallic dots 
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^ ; Abstract 

"j^ I We study the effects of the mesoscopic fluctuations on the competition between exchange and 

pairing interactions in ultrasmall metallic dots when the mean level spacing 5 is comparable or 
^ ■ larger than the BCS pairing energy A. Due to mesoscopic fluctuations, the probability to have 

o , 

CJ , a non-zero spin ground state may be non- vanishing and shows universal features related to both 

level statistics and interaction. Sample to sample fluctuations of the renormalized pairing are 
■ enlightened. 
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The competition of superconductivity and ferromagnetism and their possible coexistence 
in bulk materials has been intensively investigated in the past and it lead to a very rich 
phase diagram In addition to the suppression of superconductivity due to time reversal 
symmetric breaking, new phases with non-uniform magnetic and superconducting ordering 
were predicted 0. This picture does not hold when the size of the sample is reduced 
such that the average level spacing 6 becomes comparable with the energy scales related 
to the onset of macroscopic order. In this case fluctuation (both thermal and quantum) 
effects become very important and even the definition of ordered phase is elusive. For 
instance a small superconductor will not posses a fully developed gap ^ and the spin of a 
ferro magnet will not be macroscopically large IQ. Nevertheless distinct features reminiscent 
of the macroscopic order can be traced even in nanosized grains. 

In this Letter we investigate the competition between exchange and pairing in small metal- 
lic grains. Our approach starts from the results of a recent work of Kurland et al. who 
have shown that universal features of electrons in isolated mesoscopic grains are accounted 
for by the Hamiltonian (see also Refs. 0) 

T< = Y.Y. ^A,aCa,a + EcN'-XPf-JS'. (1) 

a=l ^={T,i} 

This description is valid in the limit of very large dimensionless conductance g = Et/S, 
being Et the Thouless energy and 6 the mean level spacing. In Eq.(|l|) the first term is the 
kinetic energy while the electron-electron interaction is expressed as a sum of three terms, 
which describe respectively charging, pairing and exchange interactions. The index a spans 
a shell of fl doubly degenerate (a = ±) time reversed single particle states of energy e^, 
and Ca,cr (cjj ^ ) are the corresponding annihilation (creation) operators. The energies 
are distributed according to the Gaussian Orthogonal Ensemble (GOE), which describes the 
case of time-reversal symmetry with no spin-orbit coupling . The interaction depends only 
on the collective variables N = c}^ a'^a^a (the number operator), T = Yl^=i Ca,-Ca,+ (pair 
creation operator), and S = Yla=i^aa^<ry(^a,a' (the total spin operator, a are the Pauli 
matrices). For isolated grains is fixed and the charging term can be ignored. Pairing 
interaction tends to favour the formation of spin singlets, while exchange tends to favour 
maximal spin, so they compete in determining the spin ordering. 

A clear picture is already emerging in the case where either pairing or exchange is present. 
In the absence of exchange interactions (J = 0), the Hamiltonian in Eq.(|l|) describes pairing 
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correlations in small metallic grains 0. Although no phase transition occurs, signatures 
of pairing correlations have been detected even in nanosized grains. For this model it 
has been shown p!0| , pJ]] that the low-energy properties are universal functions of the ratio 
6/A = 2sinh(5/A)/n (A is the BCS gap value). Upon increasing the size of the grains, 
hence decreasing the ratio 6/ A, there is a crossover [|n|] between the case of ultrasmall 
grains {6 ^ A) where pairing produces strong quantum fluctuations |TD[, and a regime 



(5 ^ A), where the BCS mean-field description remains valid |T^. Signatures of pairing 



correlations may be detectable in thermodynamic quantities [|T3[, even for ultrasmall grains. 
In absence of paring (A = 0) the properties of the model of Eq.([l|) are determined by the 
interplay between the kinetic term, which favours Pauli filling of the levels and zero total 
spin in ground state, and exchange one which tends to maximize S and eventually leads to 
the Stoner instability for J > 6. However in mesoscopic samples P, 14, 15, it is possible 
to find individual grains with a cluster of 25" closely spaced levels around the Fermi energy, 
whose ground state may have spin S even for J ^ 5. The probability Ps{J/S) of spin-S* 
ground state directly reflects the universal properties of the level statistics [0. 

In order to describe the interplay between superconductivity and ferromagnetism in small 
grains we study how the tendency to magnetic ordering is reduced when pairing is gradually 
increased. Hence, in the same spirit as in Ref. [|T7|, we consider the probability Ps of finding 
a spin-S* ground state in the regime J, X <^ 6, which is in turn related to the spontaneous 
magnetization of an ensemble of grains. All the results that will be presented are obtained 
in the half-filling scheme, where the numbers of electrons is equal to the number of levels 
Q. For an ensemble of normal grains Ps{J/S) for J ^ 5 is non-zero due to grains with 



a cluster of 25* close levels around the Fermi energy ll^. The same picture applies when 
pairing interaction is present except that, since the energy balance between spin S states 
involves the pairing correlation energy, one should take into account contributions coming 
from the entire shell of Q levels. In the weak coupling limit J,X<^6 there is a simple way to 
circumvent this problem, namely to consider a shell of 25* levels with a renormalized pairing 



coupling A2S' (see Ref. [0). Notice that whereas no sample to sample fluctuations affect 
the bare A, the actual level distribution may produce fluctuations in A2S, which we ignore 
at the moment. With this hypothesis the probability to have spin-S* in the ground state, for 
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< J - JS < A < 5, is 



where as = {S + 1){2S — l)/2, is the threshold value of the exchange below which 
the spin probability vanishes and depends on the spin and the pairing (Jj* = A2, = 
2A3/3) and the Cs are dimensionless constant depending only on the spin S (Ci = 7r^/3, 
C3/2 = Ovr^/SO) which directly reflect the universal statistical properties of the GOE level 
distribution. Compared to the case where pairing is absent the exponent as is halved. 

We now check these results with numerical calculations. We have considered grains with 
size up to = 30, 31, which are large enough to show the universal behavior of the pairing 



interaction |TT[, as in grains with much larger f2. We considered ensembles of grains whose 
single particle spectra {ea,a = realize the GOE level statistics. To this end 

we diagonalized a set of 5^x5^2 real orthogonal random matrix taking out the central fl 
eigenvalues, in order to avoid edge effects. Then we found the many-particle energies of 
Eq.(|l]) by using the Richardson exact solution [IS] for larger systems {Q > 20), whereas for 



smaller systems we used the standard numerical diagonalization of the Hamiltonian or the 
Lanczos method. For systems with even (odd) we studied the probability Pi (-P3/2) of a 
ground state with spin 1 (3/2) using a set of 10^ (10^) level configurations (systems with 
odd require more statistics because P3/2 is smaller). Results are shown in Fig.|I] and Fig.|^ 
for values of J/ 6 such that the probability that the ground state has larger spin (S* > 1 or 
S* > 3/2) is negligible. The result given in Eq.(H) reproduces quite well the numerical data, 
except for smaller 6/ A (large grains) and near J|, a regime which we discuss later. Thus 
we conclude that Pi and P3/2 show the universal features of the level statistics predicted in 
Eq.(^, namely the coefficient Cs and the power law as a function of (J — Jg). In addition 
we find a new manifestation of the universal behavior in the quantity Jg. According to the 
simple theory which leads to Eq.(^ for grains with even (odd) A^ this quantity is related 
to the renormalized two (three) levels pairing constant A2 (A3). By fitting the linear part of 
the curves in Fig.|I] and Fig.^, we find that A2 and A3 are given by a universal functions of 
6/ A, of the form A25 = 6/ \n{asS/ A) |jTO|, This result, shown in the insets of the Fig.|l| 



and Fig.0, is valid over several decades of values of the parameters. 

We now discuss more carefully the behavior near for grains with even A^ (Fig.^. 
For samples with larger pairing interaction (smaller 6 /A) the probability Pi shows a non 




J/5 

FIG. 1: The Pi{J/5) probabilities in the case $7 = 30 for the different values of the ratio 5/A 
(circles). The lines represent the fits to the data using the expression given in Eq.(|2|). In the inset, 
the numerical data of the renormalized pairing constant (circles) extracted from J* are plotted as 
a function of 5 /IS.. They coincide with the value of the renormalized pairing A2 = 5/ ln(ai(^/A) 
(with ai = 1.721). 




FIG. 2: The probabilities in the case = 31 for the different values of the ratio 5/ A. For 
convenience we plotted P^/^iJl^) = [1350P3/2(J/5)/vr4](2/^V(9^/'^) = {J - 2\-i/'i)/5. The circles 
represent the numerical data, and the lines the linear fits. In the inset, the numerical data threshold 
value as a function of the ratio 5/ IS.. The scaling behaviour is obtained as in the previous figure 
with 03/2 = 1.679. 
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FIG. 3: A typical configuration of the single-particle energy levels (we chose as an example il. = 6 
and 6/A = 25) which leads to a ground state with spin one in the tail region (left side, J/6 = 0.233) 
and in the linear region (right side, J/6 = 0.300). The Fermi level has been set to zero. 

vanishing tail for J ^ J*- We argue that this effect is due to sample to sample fluctuations 
of the renormalized A25 which we ignored in Eq.(^). To understand this point consider 
even N and samples for which the two levels in the central cluster at the Fermi energy 
are closely spaced (s2 ^ S). These samples may contribute to Pi{J/5) for J ~ Jj". The 
two-level renormalized coupling in one of these grains is determined by the configuration of 
the other levels in the VL shell and in particular depends strongly on the spacing S4 between 
next neighboring pair of levels (above and below the Fermi energy, see the right side of 
Fig.^). For most of the samples S4 ~ 35 and they will have approximately the same two- 
level renormalized coupling, which is identified as the threshold J^ in Eq.(0). However for 
a small fraction of samples we may have S4 ^ 35 (see the left side of Fig.|^) which leads 
to a smaller value of the renormalized coupling. As a consequence pairing correlations will 
be weaker and the ground state will have 5 = 1 for a value of J smaller than J^. These 
samples determine the appearance of the tail for J ^ Jl- 

To check this argument we have first analyzed the typical level configurations of samples 
contributing to Pi{J/6) (Fig.^. Indeed samples which contribute to Pi{J/6) in the tail 
region have level configurations as in Fig.|^, with a central cluster of two levels around the 
Fermi energy, and all the others very far away, whereas in the region described by Eq.(0) the 
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FIG. 4: The behavior of the spin-1 ground state probabihties for different values of the bandwidth 
0, (curves with different symbols) at the fixed value 6/ A = 25 in the tail region. The form 
Pi{J/5) = (81/647r^(5^) J(J — A)erfc[6/(J — A)] (erfc(a;) is the complementary error function) was 
chosen from the solution with three levels. The fit to the curve with = 30 (circles, in the inset) 
is obtained for X/6 ~ 1^/5 = 0.201 and b = 6.821 • lO^^^. 

other levels are closer to the two central ones. Further insight can be gained by including 
in Eq. (|^) the effect of fluctuations of a neighboring level. We consider a shell of three levels 
with two electrons. Two levels are close while the third lies far away. In this system we 
determine the approximate form of the probability distribution Pi{J/6) (caption of Fig.^, 
which is non vanishing in the tail region. Moreover we fitted this result with numerical data 
for larger systems (up to i7 = 30, see inset of Fig.p, using the value of the renormalized 
coupling as a fitting parameter (see the inset of the Fig.H). The agreement with the numerical 
data suggests that fluctuations of the renormalized pairing constant due to the statistics of 
far levels are negligible. 

We also studied the behavior of Pi{J/6) for J < J^, for fixed 5/A = 25 and for various 
values of the bandwidth Q. We used sets of GOE levels with 10'' realizations for f2 = 4, 6, 10 
and with 10^ realizations for Q = 20, 30. We see that for the values of J we could study 
Pi{J/6) becomes independent on Q for Q > 20, showing that still Pi{J/6) depends on 
features of the universal level statistics. However by decreasing J the behavior of Pi{J/6) 
will depend more and more on details of the full Q shell so universality is lost (for instance 



7 



Pi{J/S) is expected to be nonzero for J > A, the bare couphng). On the contrary by 
increasing J we recover the behaviour of Eq.(0) and all the curves collapse on the same 
line. The above analysis could be in principle carried out also in samples with odd A^. 
However the unpaired electron in the spin 1/2 ground state weakens pairing correlations. 
As a consequence the tail in P^i2{J/5) is tiny and it would require a much larger statistics 
to be investigated. 

Finally we consider grains with larger pairing interaction such that 5 < A. This is the 
crossover region to BCS-like superconductivity 0, |l^ and an analytical approach to the 
problem is a formidable task, so we studied this regime numerically. Some qualitative results 
can be inferred by looking at the evolution of the behavior of the probability distributions 
Ps{J/S) from the perturbative region 5 ^ A to the region 5 ~ A. In Fig.|^ the probability 
distributions of the smallest and largest spin ground states are shown for the systems with 
an even {Q = 30) and odd {Q = 29) number of electrons in 10^ GOE realizations. For 
smaller J we have Po{J/S) = 1 whereas for large J the maximal spin is favoured. In the 
intermediate region the stable value of S is finite but the physics is sensitive to mesoscopic 
fluctuations. Results show that this crossover region becomes narrower as the ratio 6/ A 
decreases. This trend is confirmed by the fact that, on increasing the exchange, the system 
tends to maximize the spin in the ground state, whereas in the dot regime 6 > A the 
crossover occurs gradually across all the increasing values of the spin. In the samples with 
an odd number of electrons the width of the J range between the probability of the smallest 
and largest spin probability is larger than in the even systems since the pairing is weakened 
by the presence of the unpaired electron. 

In conclusion we studied a model for the competition between superconducting and fer- 
romagnetic ordering in small metal grains. The presence of pairing implies that all the 
levels in the shell are responsible for the determination of the spin in the ground state. As 
compared to the normal case there are three new features: i) a soft threshold appears in the 
Ps directly related to the renormalized pairing coupling; ii) the power law behavior has a 
new exponent as compared to the case where pairing is absent; iii) at low enough exchange 
couplings Ps{J/S) is determined by sample to sample fiuctuations in the pairing coupling. 
By increasing both the pairing interaction and the exchange coupling we found that the 
region of the phase diagram where mesoscopic effects are important becomes progressively 
narrower. 



8 




1.25 



FIG. 5: The behavior of the probabihty distributions in the systems with an even (upper graph, 
0, = 30) and odd (lower graph, 0, = 29) number of electron, of the smallest (5 = and S = 1/2, 
decreasing curves on the left side) and largest {S = 15 and S = 29/2, curve on right side of the 
figure) spin ground states for different values of the 5/ A ratio (denoted by the numbers near each 
curve, in the case of largest spin value all the curves collapse on a single one) in a set of 10^ 
realizations. 
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